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The Viscosity of Ice. 

By B. M. Deeley, F.G.S. 

(Communicated by Henry Woodward, LL.D., F.R.S., F.G.S. Received May 18,— 

Read June 4, 1908.) 

Although the character of the motion of glaciers is now well known, and 
the velocity of their motion has, in some cases, been measured, calculations do 
not appear to have been made of the viscosity of either glacier or crystalline 
ice ; nor do there appear to be any published figures from which the viscosity 
can be quite accurately ascertained. It is possible, however, by estimating 
certain dimensions, to obtain results which, although not strictly accurate, 
give figures of the proper order of magnitude. 

The property termed viscosity is here used as defined by Maxwell.* A 
tallow candle is much softer than a stick of sealing-wax; but, if the candle 
and the stick of sealing-wax are laid horizontally between two supports, the 
sealing-wax will,’ in a few weeks in summer, bend with its own weight, while 
the candle remains straight. The candle is, therefore, a soft solid, and the 
sealing-wax a very viscous liquid. A viscous body is considered as one which 
undergoes permanent and continuous change of form under the action of a 
stress, however small. A plastic body is one which requires a definite and 
often increasing stress to produce continuous change of form. For all 
practical purposes, a substance which undergoes continuous change of form 
under very small stresses and yields at a rate proportional to the stress may 
be regarded as viscous. I have previously given some formuleef for calculating 
the viscosity of liquids, but did not then apply them to any concrete case. 
This I now propose attempting. 

If a bed of a glacier were fairly symmetrical and the ice obeyed the laws of 
viscous flow, then the ice would move faster in the middle and more slowly at 
the sides and bottom. The rate of flow was carefully measured by Tyndall, 
who placed a straight line of stakes across a glacier and found that in the 
course of a few hours the straight line became bent into a curve. We have 
no measurements which show that this curve is strictly parabolic in the case 
of any particular glacier, nor can we expect to find such, for the valleys down 
which glaciers move are generally sinuous and have irregularities in their 
beds. 

Above Trelaporte, on the Mer de Glace, crevasses form right across the ice, 
and into these crevasses the rock debris on the glacier falls. When, by 

* ‘Theory of Heat/ 1894 edition, p. 303. 
t ‘ Geol. Mag./ September, 1895, p. 408. 


The Royal Society is collaborating with JSTOR to digitize, preserve, and extend access to 
Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character. 

www.jstor.org 














Downloaded from rspa.royalsocietypublishing.org 


o 
z 
oc 


2<Z}U 


^u 
Sq.o3i/> 



Ss 


w u 

Ko 

H m 



_r a 
< z 

y 5 

b m 
2<zm 

^S=!v5^ 

h^zS 

5>ui- 

2i v u 
S0.08 m 



Ss 


w U 

Ko 

H ^ 



The Viscosity of Ice . 


251 


melting of the ice surface, these dirt bands become visible lower down, the 
dirt bands are seen to stretch across the glacier in great parabolic curves. 

The downward movement of a glacier is due to its weight and viscosity. 
In the case of a fluid of small viscosity, such as water, inertia effects cause 
the stream to be swifter in some places than can be due to the slope there, 
and slower at other places. With a glacier, the viscosity of which is enormous, 
the great viscosity enables it to transmit thrust to some considerable distance, 
and for this reason we get velocities at various points which are not wholly 
due to the effects of gravity at such points. The surface melting of glaciers 
also causes them to undergo longitudinal compression in many cases, and this 
again complicates the question. 

Tyndall measured the rate of flow of a number of glaciers, the slopes of the 
upper surfaces of which can be approximately obtained. Their thicknesses, 
however, can only be estimated. This I have done by comparing the valleys 
down which the glaciers move with other valleys of somewhat similar width 
and length. I consider it probable that the slip at the bottom is about 
15 per cent, of the surface velocity, so that the effective difference of velocity 
between top and bottom is 75 per cent, of the surface velocity. The rock 
being dense and hard, the probability of the glacier being frozen to its bed 
does not seem great. The breadth of the glaciers being much in excess of 
their thicknesses, and it being probable that the slip is much greater at the 
sides than at the bottom, the assumption has been made in the calculations 
that the glaciers are of infinite width. 

In Table I are given both the calculated viscosities for a number of glaciers 
and the figures upon which the calculations are based. These figures will be 
subjected to criticism and modified if more accurate estimates of the data can 
be made. 

As an instance of the method adopted, the calculation for obtaining 
viscosity of the Great Aletsch glacier is given below :— 

P = bodily force producing motion, being downward component of gravity 
= gravity x density x gradient 

= 980*5 x 0*92 x 65/1000 = 58*634 dynes per cubic centimetre. 


By observation, the maximum velocity V is 0*00047 cm. per sec. and the 
depth v is at the bottom of the glacier 45,000 cm. 

The equation of flow, = leads to 


fi^ = IV, and 
dr 


_ Pr 3 _ 58-634 x 45,000 2 
M 2Y 2 x 0-00047 


126,224,200,000,000, 


say 126 x 10 12 , the viscosity of water being only 0*01028 at 20° C. 
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252 Mr. R. M. Deeley. [May 18, 

The maximum shearing stress at the bottom of the glacier tending to drag 
the rocky surface along with it is the total downward component of gravity, 

Tr = 3*731 x 1476*4 = 5510 lbs. per sq. ft. 

In many cases glaciers have in my opinion* been frozen to their beds, and 
have dragged along with them, arid crumpled and contorted, the rocks over 
which they have passed. 

The shear at the bottom of the Great Aletsch glacier, amounting to nearly 
2J tons per square foot, is the average near the centre. This figure may be 
greatly exceeded at points where a rock projects above the general surface; 
for owing to the very great viscosity of the ice, thrust can be transmitted very 
considerable distances. 

Although glacier grains are capable of being deformed easily by a stress 
producing shear at right angles to the optic axis, yet from the fact that these 
grains have their axes at all angles, the motion of the glacier must be partly 
due to some other assisting cause. This I have suggested is the liquefaction 
and regelationf which takes place when molecules pass from crystal to 
crystal.^ An ice crystal may evaporate more freely on one face than 
another, thus causing a permanent and continuous transfer of molecules 
from crystal to crystal. McConnel states that a bar of ice with which he 
experimented decreased in width when the optic axis was vertical more 
rapidly than it decreased in depth. 

The effective viscosities of glaciers referred to in Table I range according 
to these calculations from 292*2 x 10 12 to 3*27 x 10 12 . It is admitted that 
the data are somewhat uncertain; but as far as rate of movement, slope, and 
thickness are concerned, the errors can scarcely be as great as will account 
for the whole of this difference. The most serious uncertainties are the extent 
to which thrust and the amount of slip at the boundaries may affect the 
results. In the case of the Great Aletsch the regularity of the slope and the 
great length of the ice stream seem to show that there is no rock basin of 
any depth below the stream, and the value of the viscosity obtained for it is 
probably the most reliable. Both the Mer de Glace and the Lower 
Grindelwald glaciers are in compression and the smaller values may be 
partly due to this, ie. to thrust giving a more rapid motion than is due to the 
slope only. In these two ice streams the glacier grains may also be smaller, 
owing to the nearness of the n4ve and the presence of ice falls, and this may 
cause the viscosity to be small. 

The Morteratsch results are also likely to be affected by thrust, especially 

* £ Geol. Mag./ 1897, p. 394. 

+ ‘ Phil. Mag.,’ February, 1888. 

% ‘ Geol. Mag.,’ May, 1895, p. 161 ; ‘ Phil. Mag.,’ May, 1895. 
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towards the end, and a mean of the three viscosities obtained for this glacier 
may he considered, to be fairly free from error. 
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Table I.—Viscosity of Glacier Ice. 


Glacier. 

Thickness, 

metres. 

Velocity. 

Slope. 

Pressure 
in dynes 
per cm. 2 

Viscosity 

Ab in. 

c.g.s. units. 

Including 

slip. 

Without 

slip. 

Inches 
per day. 

Inches 
per day. 

One in 

Mer de Glace— 







At Trelaporte . 

210 

20 

17 

13 -5 

66 *82 

29 *48 x 10“ 

At les Ponts . 

210 

23 

20^ 

12 

75 T7 

27 *51 

Above Montanvert... 

210 

26 

22\ 

11 

82 -00 

27 *64 

At Montanvert . 

210 

34 

29 

10 

90 *20 

23 '33 

Below Montanvert... 

210 

33 

28 

10 

90 -20 

24 T6 

Morteratsch, Pontre- 







sina— 







High up . 

400 

14 

12 

7 

128 '90 

292 '2 

Middle. 

290 

11 

9 

10 

90 -20 

143 '4 

Above end . 

190 

7 

6 

10 

90 -20 

92 '31 

Lower Grindelwald ... 

100 

22 

18| 

25 

36 *09 

3-274 

G-reat Aletsch. 

450 

— 

16 

15 ‘4 

58 '63 

126 -2 
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A viscosity of 125 x 10 12 appears to be as near an estimate as can yet be 
made of the average effective viscosity of a Swiss glacier. In winter the 
viscosity is probably double this figure owing to the lower temperature. 

J. C. McConnel* appears to have been the first to observe the exact condi¬ 
tions under which ice is capable of being deformed without fracture by stress. 
He showed that a crystal of ice can be sheared by very small stresses in a 
direction at right angles to the optic axis, and that the rate of shear becomes 
greater as the stress is increased. His two main conclusions are (1) that the 
friction between the particles of ice along the shear planes becomes greater as 
the temperature falls, (2) that when the molecules of ice slide over each 
other the cube of the friction varies as the square of the velocity. It must 
be remembered that McConnel died before his paper was quite finished, and 
that, therefore, his calculations were not completed, nor was his paper 
subjected to final correction. It will, consequently, not be out of place to 
reconsider his experimental results in detail. 

In Table II, I have set out a number of the experimental results he 

* ‘ Boy. Soc. Proc., J May, 1891, p. 323. 

VOL. LXXXI,—A. S 
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254 Mr, R. M. Deeley. [May 18, 

obtained with bars of ice loaded with various weights. Although in one place 
only he states the distance between the supports of the ice bars, it is clear 
that the distance was 51 mm. He used a solid iron frame placed in a small 
box to keep the temperature from varying rapidly, and to check evaporation. 
In no case is the full length of the ice bar given; this I have assumed to be 
60 mm. He states that during the experiments the bars lost weight by 
evaporation, and as he gives the rate of evaporation, I have calculated the 
probable size of the bar at each reading of deflection. This can readily be 
done, as McConnel gives the size of the bar before and at the end of the 
experiment, and also gives the times when the weights were altered and the 
deflections were measured. In a few cases the figures he gives for the dura¬ 
tion of the tests are slightly in error; these have been corrected. 


Table II.—Experiments made by J. C. McConnel.* 


Time, 

hours. 

Bending, 

mm. 

| Weight, 

| kilos. 

t 

Breadth, 

mm. 

Depth, 

mm. 

Temp., 

°C. 

Shear F., 
dynes per cm. 2 

Viscosity, 
c.g.s. units. 

0 -5000 

0*118 

2*5 

14 *61 

32*27 

- 6*7 

68 *5 x 10 4 

22 *64 x 10 1 ® 

0 -5167 

0*089 

2*5 

14 *59 

12 *26 

- 7*2 

68 *6 

31 *07 

1*167 

0*086 

1 *445 

14 *57 

12 *24 

- 6*7 

39 *8 

42 *17 

1*167 

0*128 

1 *445 

14 *53 

12*22 

- 6*7 

40 *0 

28 *45 

1*867 

0T25 

0 *915 

14 *44 

12 *16 

- 8*9 

25 *7 

29 *90 

1 *483 

0*086 

0*915 

14 *37 

12 *12 

-10 *0 

25 *9 

34 *81 

14 *83 

0*086 

0*295 

14 *30 

11 *93 

-15 *0 

8*7 

117 *20 

0*8 

0*328 

2*5 

13 *84 

11 *76 

- 8*9 

75 *4 

14 *35 

0*7 

0*032 

2*5 

11 *56 

13 *50 

-15 *3 

78 *6 

134 *20 

2*333 

0*272 

2*5 

11 *51 

13 *47 

-11 *1 

79 *2 

52 *99 

0 *6833 

0*260 

2*5 

11 *45 

13 *44 

- 7*0 

79 *7 

16 *35 

0 *2667 

0*157 

2*5 

11 *44 

13 *43 

- 5*6 

79 *9 

10 *59 

0*275 

0*257 

2*5 . 

11 *43 

13 *42 

- 3*9 

80 *0 

6*68 

0 *1916 

0*257 

2*5 

11 *42 

13 *41 

- 2*8 

80 *1 

4*66 

0 *1083 

0*164 

2*5 

11 *42 

13*41 

- 2*2 

80 *1 

4 *13 

0 *1083 

0*204 

2 *5 

11 *41 

13 *41 

— 1 *7 

80 *2 

3 *32 

0 *2416 

0*310 

2*5 

11 *38 

13 *39 

- 61 

80 ‘5 

4*90 

0*4083 

0*304 

2*5 

11 *23 

13 *29 

-13*0 

82 *2 

8*62 

0 *6416 

0*440 

2*5 

11 *21 

13 *27 

-13 *2 

82 *5 

9*39 

3*917 

0*146 

2*5 

16 *03 

17 *73 

- 5*0 

43 *2 

90 *50 

1*408 

0 *227 

2*5 

15 *95 

17 *67 

- 5*0 

43 *6 

21 *10 

1*542 

0*084 

1 *29 

15 *90 

17 *63 

- 6*7 

22 *7 

32 *47 

0 *4583 

0 073 

2*5 

15 *86 

17 *61 

- 7*5 

44*0 

21 *55 

0*4 

0*119 

3*79 

15 *84 

17 *60 

, - 8 *0 

66 *8 

17 *51 

1 *058 

0*342 

3*79 

15 *82 

17 *58 

- 8*9 

66 *9 

16 *15 

0 ‘7416 

0*146 

2 *5 

15 *79 

17 *56 

-10 *0 

44*3 

17 *56 

12 *83 

0*288 

0*62 

15 *56 

17 *41 

- 7*8 

11 *3 

39 *39 

0 *5916 

0T65 

2*5 

15*33 

17 *26 

- 6*1 

46 *4 

12 *99 

0 *4083 

0*235 

3 *79 

15 *31 

17 *24 

- 5*9 

70 *5 

9*56 


# See ‘ Roy. Soc. Proc.,’ March 12, 1891, p. 323. 


Fig. 1 shows a bent bar of ice resting upon its supports with the optic 
axis vertical. When the bar is bent the ice assumes the form given by the 
dotted lines at the point of flexure ; but for all practical purposes it may be 
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assumed that the bending is as shown by the full lines. This assumption 
much simplifies the calculations, and does not appear to me to introduce any 
serious error. 



f'*- F/c. 2. 


When, as shown in fig. 2, a mass of a viscous substance is distorted by the 
application of a stress to the plane XX, and the plane YY is fixed, the 
distortion is as shown by the dotted lines, and the stress at all planes between 
and parallel to the planes XX and YY is the same. Now, according to fig. 1, 
the ice bar is distorted by shear in exactly the same way, and I have regarded 
the weight on the centre of the bar as producing a stress parallel with and 
coincident with the upper and lower faces of the bar. 


Calculations for Viscosity . 

I = distance between supports, ‘cen- A = total area of horizontal sur- 
timetres. face, square centimetres. 

8 = deflection, centimetres. h = time, hours. 

r = thickness, centimetres. W = load, dynes. 

The shear cr is 8/^lh per hour ; the total shearing force in the cross-section 
is yuA'cr/60 2 , where A' is the area of the cross-section. This must be equal 
to JW, so that 

900W Ih 


P : 


A'8 


For the ninth test of Table II, 

l = 5T ; 8 = 0*0032 ; r = 1*35 ; A' = 1*35 x 1*156; k = 0*7; 
W = 2500 x 981 = 2,452,500, 
so that p = 134*2 x 10 10 , 

W _ 2,452,500 _ 


and the shearing force = 


2A y 2 x 1*561 


= 78*6 x 10 4 . 


In the last two columns of Table II are given the shearing stresses to 
which the bars were subjected and also the viscosities. 

s 2 
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256 Mr. R M, Deeley. [May 18, 

On Diagram I these viscosities are plotted, and against each plotting 
is marked the shearing stress (in dynes per square centimetre x 10 4 ) to which 
the ice was subjected. 

The equation to the curve on the diagram is 

login ji = 0*301 + 0*153£—0*0023l£ 2 , 

so that log /6 is parabolic, where fju is the viscosity in dynes per square 
centimetre x 10 10 and t is the temperature below zero C (considered positive). 



Diagram I. 


gglt will be noticed that out of 29 tests, one of which is stated by 
McConnel to be inaccurate, and two more of which also appear unreliable, 
leaving 26 presumably accurate results, 14, or more than half, give results 
agreeing within 5 units with the curve. The want of exact agreement with 
the curve seems to be mainly due to the varying temperatures experienced 
during the continuance of the experiments. 

McConneTs experiments show that when the weight was taken off the ice- 
bar, in addition to purely elastic rise it recovered slowly some of the bending 
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1908.] The Viscosity of Ice . 257 

the weight had produced. Many, if not all, soft solids do this. The recovery 
was very slow, and almost ceased after a time. 

It will be seen, on referring to Diagram I, that the experiments made with 
the small stresses give, in some slight degree, the greatest viscosities, but 
this is not true to anything like the extent McConnel thought. In this 
connection he states that when the weight is changed the alteration in the 
rate of depression is great out of all proportion, e.g ., the alteration in the 
rate of bending from 0*0058 to 0*41 when the weight is changed from 
0*295 to 2*5 kilos. In this case, however, there was a change of tem¬ 
perature from —15° C. to — 8°*9, and, as'will be seen from the diagram, the 
change of temperature will account for almost the whole of the change of 
viscosity. He should have compared the case in which the temperature was 
— 15° with the one in which it was —15°*3 C., and he would have found that, 
with stress varying in the proportion of 1 to 9, the rate of shear was very 
nearly proportional to the stress. 

Mugg found that continuous shear could not be produced in some samples, 
of ice when very lightly stressed. 

It may be that the small departure from the rate of shear being pro¬ 
portional to the stress, when the stress is small, is to some extent due to 
want of regularity in the structure of the ice tested, and may al$o be in some 
measure connected with the property ice possesses of slightly recovering its 
original form after distortion. 

There being some uncertainty as to the actual behaviour of other highly 
viscous substances under stress, I made some experiments with a bar of pitch. 
It was shaped by heating and pressing into a bar 150 mm. long and 14 mm. 
square, and placed upon supports 117 mm. apart. The load at the centre of 
the bar was 98*5 grammes, partly due to the weight of the bar itself and the 
finger used to indicate the bending. After the application of the weight, the 
bar was allowed to bend until the finger reached the scale. The readings 
then taken are given in Table III. At 12.19 p.m. a weight of 67*6 grammes 
was removed from the centre of the bar, which was then stressed onty by the 
finger resting upon it and its own weight. After the removal of the weight, 
readings were taken for a time every 15 seconds. On Diagx^am II some of 
the readings taken before and after the weight was removed are plotted. 
The diagram shows that there was an elastic unbending of the bar, followed 
by a further slow unbending, which ceased in about 4 minutes. The bar 
then again began to yield viscously under its own weight. 

From this it is concluded that a viscous substance under stress requires 
time after a change in the magnitude of the stress has been made for the 
rate of shear to become proportional to the altered stress. Also, that a 
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258 The Viscosity of Ice . 

viscous (liquid) substance may be elastic, brittle, and may be able to partly 
regain its shape after being deformed by viscous distortion. 


Table III.—Sample of Pitch. April 17, 1908. Bar 150 mm. long, 
14 mm. broad, 14 mm. deep. Supports 117 mm. apart. 


Time. 

Deflection, mm. 

Weight. 

Temperature, ° C. 

11.30 

0*50 

On 

12 *8 

11.47 

1 -oo 

77 

12 *8 

12.1 

1 *37 

77 

12 *9 

12.10^ 

1 *62 


12 *9 

12.19 

1*87 

77 

12 *95 

12.19i 

1*77 

Off 

12 *9,5 

12.20 

1-75 

97 

12 *95 

12.20£ 

1*73 

77 

12 *95 

12.20f 

1*71 

77 

12 *95 

12.23 

1*69 

77 

12 *95 

12.31 

1 *70 

77 

13*0 

12.56 

1 -71 

77 

12 -95 

1.53 

1 -75 


12 *75 

3.11 

1*92 

77 

12 *8 


During the experiment with pitch, the thermometer and experimental 
apparatus were protected by a glass case and duster. 



Diagram II. 


Since writing the above, my attention has been called to the experiments 
made by Trouton* on several viscous substances, including pitch. He states 
that “it was noticed that there was a slow partial movement towards 
recovery on removal of the force of traction, which gradually fell to zero with 
time.” 

We may, I think, fairly conclude, from a consideration of McConnel’s 
experiments and these above described on pitch, that the nature of the shear 


* ‘ Boy. Soc. Proe./ vol. 77, p. 429. 
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which can be produced in an ice-crystal at right angles to the optic axis very 
closely obeys the laws of viscous flow. If this view should be sustained, it is 
not only true, as Professor H. A. Miers shows, that a liquid of low viscosity 
may have a crystalline structure, but also that, as in the case of crystalline 
ice, a solid may be liquid along one plane only. 

Prom the equation to the curve on Diagram I, the viscosity of an ice- 
crystal in a direction at right angles to the optic axis is about 2 x 10 10 at the 
freezing point. At this temperature the viscosity of a glacier is about 
12,500 x 10 10 , or 6250 times as great. The effective viscosity of a glacier is 
therefore due in a great measure to some other consideration. The optic axes 
of glacier grains are at all angles and they lock each other. 

I am much indebted to Mr. P. H. Parr for the assistance he has rendered in 
calculating the results from the figures given. 


Vortices in Oscillating Liquid . 

By Lord Rayleigh, O.M., Pres. R.S. 

(Received July 25, 1908.) 


In a paper “ on the Circulation of Air observed in Kundt’s Tubes, and on some 
Allied Acoustical Problems,”* I applied the equations of viscous incompressible 
fluid to show that the effect of the bottom of the containing vessel was to 
generate permanent vortices in the vibrating fluid. It was remarkable that 
the intensity of the vortical motion, when fully established, proved to be 
independent of the magnitude of the viscosity, so that the effects could not 
be eliminated by merely supposing the viscosity to become extremely small. 
The expression found for the vortices was simple. The horizontal component 
u of the primary motion near the bottom being u = cos kx cos nt , the 
component velocities of the vortical motion are 


u = 


, Uq 2 sin 2 kx 


-2ky 


(1-2 ley), 


, u 0 2 cos 2 kx 

v— 4 


~ 2k y 2ky, 


y being measured upwards from the bottom, and V ( = n/k) the velocity of 
propagation of waves of the length in question. According to these 
* ‘Phil. Trans.,’ vol. 175, p. 1, 1883 ; ‘Scientific Papers,’ vol. 2, p. 239. 











